WFS Maps

Setup

library(dplyr)
library(sf)
library(readr)
library(mapview)
mapviewOptions(fgb=FALSE)
library(Rspatialworkshop)
library(ggplot2)
library(leaflet)
library(httr)

Queyring WFS

Just to show scope of different things we can do linking to REST Services as well - here are links to a couple examples- spatial wfs services and accessing arcgis rest services

We can actually use REST services for equivalent WMS layers we displayed in maps above! See National Map Services page

We can get a listing of all ESRI REST services here

Query all Native Lands WFS

url <- parse_url("https://services.arcgis.com/P3ePLMYs2RVChkJx/arcgis/rest/services")
url$path <- paste(url$path, "Native_Lands_New/FeatureServer/2/query", sep = "/")
url$query <- list(where = "1=1",
                  outFields = "*",
                  returnGeometry = "true",
                  f = "geojson")
request <- build_url(url)


NativeLand <- read_sf(request)

mapview(NativeLand, col.regions='green',map.types = "OpenStreetMap")

Query Trails WFS

REST service and a query for ‘Pacific Crest National Scenic Trail’ - bonus add a default topography brackground with it

url <- parse_url("https://carto.nationalmap.gov/arcgis/rest/services")
url$path <- paste(url$path, "transportation/MapServer/37/query", sep = "/")
url$query <- list(where = "Name = 'Pacific Crest National Scenic Trail'",
                  outFields = "*",
                  returnGeometry = "true",
                  f = "json")
request <- build_url(url)

PCT <- read_sf(request)

mapview(PCT, color='dark green',map.types = "OpenTopoMap")